Bottomonia suppression in heavy-ion collisions 
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Using the two-component model that includes both initial production from nucleon-nucleon hard 
scattering and regeneration from produced quark-gluon plasma (QGP), we study the effect of 
medium modifications of the binding energies and radii of bottomonia on their production in heavy- 
ion collisions. We find that the contribution to bottomonia production from regeneration is small 
and the inclusion of medium effects is generally helpful for understanding the observed suppression 
of bottomonia production in experiments carried out at both the Relativistic Heavy Ion Collider 
(RHIC) and the Large Hadron Collider (LHC). 
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I. INTRODUCTION 



Since the suggestion by Matsui and Satz [l| that sup- 
pressed production of J in relativistic heavy-ion colli- 
sions could be a signature for produced QGP, there have 
been many experimental [l, Q and theoretical studies [3- 

ton this very interesting phenomenon; see, e.g., Refs. 
[l3| for a recent review. Although their original idea 
was that the screening of color charges in the QGP would 
prohibit charm and anticharm quarks from forming the 
J/tp and thus suppress its production, lattice QCD calcu- 
lations of the J / ip spectral function showed, on the other 
hand, that the J/ip could survive inside QGP up to the 
so-called dissociation temperature [13, [l^- However, the 
dissociation temperature depends on how the J/'0 spec- 
tral function is extracted from lattice data and it is not 
yet clear if it is far from the critical temperature jlGj . 
Studying J/ip suppression in relativistic heavy- ion colli- 
sions can thus provide not only a signature for the QGP 
but also a probe of its properties @,H, [13, [3 • ^ quanti- 
tative study of J/ip production in heavy-ion collisions is, 
however, complicated by their absorption in the initial 
cold nuclear matter and regeneration in the QGP. Since 
the T is a more strongly bound state of bottom and anti- 
bottom quarks than the J/ip, its production in heavy- ion 
collisions is expected to be less affected by initial cold nu- 
clear matter effects. Furthermore, the much smaller num- 
ber of bottom quarks than charm quarks that is produced 
in heavy-ion collisions makes the contribution of regener- 
ation from the QGP to T production also less important. 
Therefore, studying T production in heavy-ion collisions 
would provide a cleaner probe of the properties of QGP 
and also the in-medium properties of bottomonia [l9l[20j . 
Recently, the nuclear modification factor Raa of the sum 
of bottomonia T(1S), T(2S) and T(3S), defined by their 
yields relative to those from p+p collisions multiplied by 
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the number of initial binary collisions, in Au-|-Au colli- 
sions at y/sNN = 200 GeV was measured by the STAR 
Collaboration at RHIC [IJ, while that of T(1S) [H] and 
the relative suppression of T(2S) and T(3S) to T(1S) [1^ 
in Pb-|-Pb collisions at ^/sWn = 2.76 TeV were measured 
by the CMS Collaboration at LHC. In both experiments, 
the measured Raa was seen to decrease with increas- 
ing centrality of collisions. In the present study, we use 
the two-component model [^-[Hl, which was previously 
used for studying J/ip production in heavy-ion collisions, 
to show that these experimental results allow us to ob- 
tain useful information on the properties of bottomonia 
in QGP. 

The paper is organized as follows. We first briefly re- 
view in Sec. |TT] the two-component model for bottomo- 
nia production in heavy-ion collisions. We then compare 
in Sec. Illll the calculated nuclear modification factors of 
bottomonium from the model with experimental data. 
Finally, a conclusion is given in Sec. IIVI 



II. THE TWO-COMPONENT MODEL 

In the two-component model, bottomonia are pro- 
duced from both initial nucleon-nucleon hard scattering 
and regeneration in the produced QGP. The numbers 
of initially produced bottomonia arc proportional to the 
number of binary collisions among nucleons in the two 
colliding nuclei. Whether these bottomonia can survive 
after collisions depends on many effects from both the 
initial cold nuclear matter and the final hot partonic and 
hadronic matters. For the initial cold nuclear matter 
effect, the shadowing and the nuclear absorption effect 
are included. The shadowing is the modification of the 
parton distribution in a nucleus from that in a sirigle 
nucleon. In this study, wc use the EPS09 package [2^ 
and assume that the effect is proportional to the path 
length of a parton in nucleus [III, [2^1 ■ Using the exper- 
imental data on the average transverse momentum of T 
[2^ [26j . the ratios of gluon distribution function in a 
nucleus to that in a single nucleon at a; = rriT / ^/smn 
are 1.12 and 0.92 for collisions at RHIC and LHC en- 
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ergies, respectively, with the ratio larger than one called 
the anti-shadowing. For the nuclear absorption cross sec- 
tion of bottomonium in cold nuclear matter, its value is 
about 4.3 mb at RHIC based on the recently measured 
nuclear modification factor i?dAu = 0.78 ±0.28 of midra- 
pidity bottomonium in d+Au collisions at RHIC j27l |. 
This value is larger than the value of 2.8 mb for the nu- 
clear absorption cross section of charmonium in heavy- 
ion collisions at RHIC [1^ . The nuclear absorption cross 
section of bottomonia at LHC is, however, not known em- 
pirically. Since the nuclear cross section of bottomonium 
is expected to decrease with increasing collisions energy, 
we take the value 4.3 mb as an upper limit at LHC and 
consider the two cases of and 4.3 mb. For thermal 
dissociation of bottomonia, it includes the initial dissoci- 
ation in the produced QGP of temperatures above their 
dissociation temperatures and the subsequent dissocia- 
tion by quarks and giuons in the QGP and hadrons in 
the hadronic matter (HG). 

For the local initial temperature, we estimate it 
from the local entropy density ds/dr] = C[(l — 
a)'T'part/2 -I- a ncoii] [29| in terms of the number den- 
sities of participants and binary collisions npart(coU) = 
AiVpart(coii)/(ToAa;A?/) with A7Vpai.t(coU) being the num- 
ber of participants (binary collisions) in the volume 
ToAa;Aj/ of the transverse area AxAy and tq being the 
initial thcrmalization time. For the value of tq, we use 
0.9 and 1.05 fm/c from the viscous hydrodynamics [s^l 
for Au+Au collisions at y/s^N ~ 200 GeV at RHIC and 
Pb+Pb collisions at ^/snn — 2.76 TeV at LHC, respec- 
tively. The parameters C and a are, respectively, 0.11 
and 18.7 for RHIC and 0.15 and 27.0 for LHC from fit- 
ting measured charged particle multiplicities j7l| . With 
the quasiparticle model of three flavors for the equation 
of state of QGP and the resonance gas model for that 
of HG 0, [m , the initial maximum temperatures are 324 
and 390 MeV and mean temperatures are 269 and 311 
MeV for central collisions at RHIC and LHC, respec- 
tively. After initial thcrmalization, a schematic viscous 
hydrodynamics [ill: 113 then used to describe the ex- 
pansion of the hot dense matter with the specific shear 
viscosity rj/s taken to be 0.16 and 0.2 for the QGP at 
RHIC and LHC [s^, respectively, and 0.8 for the HG 

M- 

The thermal dissociation by partons in the QGP or 
hadrons in the HG of bottomonia that have survived from 
initial dissociations is described by the rate equation for 
the number Ni of bottomonia of type i 
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where r is the longitudinal proper time and P^ is the 
thermal decay width of bottomonia. The number of 
equilibrated bottomonia of type i is given by iV°'^ = 
7^i? rii fV0{Ti — T), where rii is its number density in 
grandcanonical ensemble; / is the fraction of QGP in 
the mixed phase and is 1 in the QGP; and 6{Ti — T) 
is the step function with Ti being the dissociation tem- 



perature of bottomonia of type i. The chemical and 
kinetic off-equilibrium of bottom quarks in the QGP is 
included through their fugacity 7 and relaxation factor 
R, respectively, with the former obtained from the con- 
servation of bottom flavor [ll| and the latter defined as 
R{t) = 1 — exp[— (t — To)/req] where Teq is the relaxation 
time of bottom quarks. With Teq ~ mg [s^] and Teq =4 
fm/c for charm quarks [llj . the relaxation time for bot- 
tom quarks is about 14 fm/c. We note that the second 
term on the right hand of Eq.([T]) takes into account the 
regeneration of bottomonia from bottom and antibottom 
quarks in the QGP. 
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FIG. 1: Thermal decay widths of bottomonia as functions 
of the temperature of QGP without (upper panel) and with 
(lower panel) medium effects. 

For the thermal decay widths of bottomonia in QGP, 
we calculate them up to the next-to-leading order (NLO) 
in perturbative QCD (pQCD) [3^. While in the leading 
order (LO) a bottomonium is dissociated by absorbing a 
thermal gluon, in the NLO it is dissociated by the gluon 
emitted from a quark or gluon in the QGP. The squared 
invariant amplitudes for these processes are the same as 
those given in Ref. |ll| for charmonia except the heavy 
quark mass and can be found in Ref. [Tl1 | . In term of the 



resulting bottomonium dissociation cross section a'^^^^ , 
the thermal decay width of a bottomonium is given by 



r(T) = Y.j 



d^k 
(2^ 



t-,ei(fc)n,(fc,T)af'-(fc,T), (2) 



where i denotes the quarks and giuons in the QGP; rii is 
the number density of parton species i in grand canon- 
ical ensemble; and Vj-d is the relative velocity between 
the scattering bottomonium and parton. For the ther- 
mal width in the mixed phase, it is taken to be a lin- 
ear combination of those in the QGP and the HG, i.e., 
r(T,) = / rQGP(T,) + (1 - /)rHG(r,), where / is the 
fraction of QGP in the mixed phase. For the dissocia- 
tion cross section of T(1S) in the hadron gas (HG), we 
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use the factorization formula [35| and assume that those 
of excited bottomonia are proportional to their squared 
radii. In the upper panel of Fig. [1] we show the ther- 
mal decay widths of bottomonia in QGP calculated with 
their masses and radii in free space, i.e., 9.5 GeV and 
0.23 fm, 9.97 GeV and 0.59 fm, 10.19 GeV and 0.95 fm, 
9.9 GeV and 0.46 fm, and 10.25 GeV and 0.82 fm for 
T(15), T{2S), TiSS), Xb(lP), and Xb(2P), respectively, 
obtained from the Cornell potential with the vacuum 
screening mass ^ =0.18 GeV [11]. These thermal decay 
widths of bottomonia are appreciable and increase with 
increasing temperature. We note that the thermal de- 
cay widths of bottomonia in HG are significantly smaller 
than those in QGP. 
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FIG. 2: Binding energies (upper panel) and radii (lower panel) 
of bottomonia as functions of the temperature of QGP. 

To include medium effects on the properties of bot- 
tomonia in QGP, we consider the modification of the po- 
tential energy between bottom and antibottom quarks 
due to the Debye screening of color charges. This is 
achieved by using the screened Cornell potential [11] : 



Vir.T) 



1 



a 



(3) 



with (T = 0.192 GeV^ and a = 0.471. The screening mass 
/i(r) depends on temperature, and we use the one given 
in pQCD, i.e., ii{T) = y^Nj3 + Nf/6 gT, where iV^ and 
Nf are numbers of colors and light quark flavors, respec- 
tively. We note that compared to results from the lattice 
QCD calculations for a heavy quark and antiquark pair 
in the QGP [37|, this potential is closer to their internal 
energy at high temperature but to their free energy at low 
temperature. The screened Cornell potential thus inter- 
polates smoothly the expected temperature-dependent 
potential between heavy quark and antiquark |38| . The 
binding energies and radii of bottomonia in the QGP can 
be obtained by solving the resulting Schrodinger equa- 
tion for the bottom and anti-bottom quark pair. Taking 



their masses to be mi, = 4.746 GeV (36[ and the QCD 
coupling constant g = 1.87 as in our previous study of 
J/ip suppression in heavy- ion collisions [Tlj . the results 
are shown in the upper and lower panels of Fig. [5] It is 
seen that the dissociation temperatures of T(1S), T(2S), 
T(3S), XbilP) and Xb(2P) are 4.0, 1.67, 1.12, 1.51 and 
1.09 Tc, respectively, and their radii increase with increas- 
ing temperature. The thermal decay widths of bottomo- 
nia obtained with their in-medium binding energies and 
radii are shown in the lower panel of Fig. [TJ and they in- 
crease with temperature and diverge at their dissociation 
temperatures. Compared to those obtained with binding 
energies and radii in free space, medium effects enhance 
the widths significantly. We note that the inclusion of 
thermal decay widths of heavy quarkonia effectively takes 
into account the imaginary part of the potential energy 
between heavy quark and antiquark at finite temperature 

MM. 



III. RESULTS 

To calculate the nuclear modification factor of bot- 
tomonia in heavy-ion collisions requires information on 
their numbers produced in p+p collisions at same en- 
ergy. Since the numbers of bottomonia produced in p+p 
collisions at LHC are not available, we use in the present 
study those of T(1S), T(2S) and T(3S) from the experi- 
mental data in j>+p collisions at ^Jsnn — 1-8 TeV mea- 
sured by the CDF Collaboration at the FermiLab |4]| . 
and those of Xb(lP) and Xb(2P) from their contributions 
to T(1S) [121 based on the branching ratios of about 0.24 
and 0.13, respectively. For p-l-p collisions at RHIC, the 
numbers of bottomonia were not individually measured, 
so we use I]n=i~3 B{nS) x da/dy\y^o{nS) = 114 pb [27| . 
where B(nS) and da / dy\y^o(nS) are, respectively, the 
branching ratio and differential cross section in rapid- 
ity for T(nS), and assume their relative abundances are 
the same as those at the LHC. For the initial number of 
bottom quark pairs in determining the fugacity of bot- 
tom flavor, it is obtained from da^^ / dy\y=o = 1.34 fih 
for RHIC [i^ and cr^^ = 17.6 iih in the rapidity range 
|y| < 0.6 from the CDF Cohaboration for LHC. Also, 



we need the contribution of excited bottomonia to T(1>S') 
in p-t-p collisions at same energy. Since this has not been 
measured, we use the information obtained from p-t-p 
collisions at ^/snn = 1-97 TeV by the CDF Collabora- 
tion at the FermiLab [42| as they are known to be es- 
sentially independent of the collision energy [4l|, i.e., the 
contributions from Xb{lP), Xb{2P), T{2S) and T(3S') to 
T(ls) are taken to be 27.1, 10.5, 10.7 and 0.8 %, respec- 
tively. The resulting nuclear modification factors (-Raa) 
of bottomonia obtained without the cold nuclear effect 
in Au-I-Au collisions at -y/sjvw = 200 GeV at RHIC and 
in Pb-|-Pb collisions at ^/snn = 2.76 GeV at LHC are 
shown in Fig. [3l It is seen that the Raa of directly pro- 
duced T(1S) is close to one even in central collisions while 
those of excited bottomnoia are small. 
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FIG. 3: Nuclear modification factor -Raa of bottomonia in 
Au+Au collisions at -y/ijvjv ~ 200 GeV at RHIC (upper 
panel) and in Pb+Pb collisions at -yijvjv = 2.76 GeV at LHC 
(lower panel) without the inclusion of the cold nuclear matter 
effect 



In the upper panel of Fig. |31 the calculated Raa of 
the sum of T(1S), T(2S) and T(3S) in Au+Au collisions 
at y/SNN = 200 GeV at RHIC as a function of the par- 
ticipant number is shown and compared with the exper- 
imental data from the STAR Collaboration 211. The 



solid line is obtained without the cold nuclear matter ef- 
fect, while the upper dotted line includes only shadowing 
effect and the lower dotted line includes both the shad- 
owing and the nuclear absorption effect. It is seen that 
because of the large experimental errors, results from our 
model with and without the cold nuclear matter effect all 
can describe the data for all centralities from the RHIC. 

For Pb+Pb coUisions at y^sKw = 2.76 TeV at LHC, 
only the Raa of T(IS') has been measured by the CMS 
Collaboration Our results for the Raa of T(1S) 

without the cold nuclear matter effect is shown by the 
solid line in the lower panel of Fig. 2] Compared with 
that measured by the CMS Collaboration [22l, these re- 
sults agree with the data for both peripheral (20-100%) 
and most central (0-10%) collisions. For midcentral (10- 
20%) collisions, our model significantly underestimates 
the measured Raa- Similar results for LHC were also ob- 
tained in Ref. j45j based on the bottom and antibottom 
quark potential that was taken to be their internal en- 
ergy from the lattice QCD. Also shown in the lower panel 
of Fig. |4] are results obtained by including the shadowing 
effect (upper dotted line) and also the nuclear absorption 
with a cross section of 4.3 mb (lower dotted line). It is 
seen that the inclusion of the cold nuclear matter effect, 



FIG. 4: Nuclear modification factor Raa of the sum of T(1S), 
T(2S) and T(3S) in Au+Au collisions at ySlTiv = 200 GeV 
at RHIC (upper panel) and that of T(1S) in Pb+Pb collisions 
at ^/snn = 2.76 TeV at LHC (lower panel) as functions of 
the participant number. Solid and dashed lines are, respec- 
tively, results with and without medium effects on bottomo- 
nia. While the upper dotted line is the result including also 
the shadowing effect, the lower dotted line further includes 
the nuclear absorption with a cross section of 4.3 mb. Exper- 
imental data are taken from Refs. [2ll. [221 . 



particular the nuclear absorption, leads to too small bot- 
tomonia Raa compared with the experimental data. We 
note that most of the suppression of T(1S) comes from 
those of its excited states, as seen from the results shown 
in Fig. 131 Also, the contribution from regeneration the 
the Raa of bottomonia is less than 1% at both RHIC and 
LHC and for all centralities as a result of the small num- 
ber of bottom quarks and the much longer bottom quark 
relaxation time than the lifetime of produced QGP. For 
comparison, the i?AA's for the case without medium ef- 
fects on bottomonia are shown by dashed lines in Fig. 21 
and they are larger than those with medium effects as ex- 
pected. In this case, the calculated Raa of T(IS') is large 
compared to the experimental data from RHIC, partic- 
ular for more central collisions. This is also the case for 
heavy ion collisions at the LHC except for midcentral 
collisions where the result obtained without medium ef- 
fects can better describe the experimental data. It is not 
clear if this indicates a change of the bottomonia sup- 
pression mechanism in midcentral collisions. Improved 
experimental data are essential for resolving this puzzle. 
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IV. CONCLUSION 

In conclusion, using the two-component model (llj 
that includes both initial production from nuclcon- 
nucleon hard scattering and regeneration from produced 
quark-gluon plasma, we have studied bottomonia pro- 
duction in heavy-ion collisions at RHIC and LHC by in- 
cluding the medium effects on the thermal properties of 
bottomonia and their dissociation cross sections. With 
the expansion dynamics of produced hot dense matter 
described by a schematic viscous hydrodynamics and in- 
cluding the thermal dissociation of bottomonia as well as 
the regeneration of bottomonia by using a rate equation, 
our model describes successfully the experimental data 
from RHIC and reasonably those from LHC on bottomo- 
nia suppression. Our results indicate that the contribu- 
tion of regenerated bottomonia is small. We have also 
studied the cold nuclear matter effect due to the shadow- 
ing or anti-shadowing of the parton distribution function 
in the nucleus and the nuclear absorption of bottomo- 
nia on their Raa ■ Although the anti-shadowing effect at 
RHIC enhances the Raa, the nuclear absorption reduces 
it. Present experimental data at RHIC are, however, not 



accurate enough to differentiate these effects. For the 
Raa of bottomonia in heavy- ion collisions at LHC, the 
experimental data arc found to be better described by 
the absence of nuclear absorption, consistent with the 
expected decreasing effect of nuclear absorption with in- 
creasing collision energy. Furthermore, the present study 
shows that the inclusion of medium effects on bottomonia 
is essential for describing the experimental observations 
at RHIC aw well as at LHC except for midcentral col- 
lisions where results without medium effects can better 
describe the data. More accurate data from future exper- 
iments are needed to obtain more definitive information 
on the properties of bottomonia in the QGP. 
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